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ABSTRACT 



The acceleration of solar flare ions during magnetic reconnection is explored 
via particle- in-cell simulations that self-consistently and simultaneously follow the 
motions of both protons and a particles. We show that the dominant heating of 
thermal ions during guide fleld reconnection, the usual type in the solar corona, 
results from pickup behavior during the entry into reconnection exhausts. In 
contrast to anti-parallel reconnection, the temperature increment is dominantly 
transverse, rather than parallel, to the local magnetic fleld. A comparison of 
protons and alphas reveals a mass-to-charge {M/Q) threshold in pickup behavior 
that favors heating of high-M/Q ions over protons, which is consistent with 
impulsive flare observations. 

Subject headings: acceleration of particles — magnetic reconnection — Sun: 
corona — Sun: flares 



INTRODUCTION 



The generation of energetic particles during flares remains a central unsolved issue in 
solar physics. Extensive observational evidence indicates that a substantial fraction of the 
energy released during a flare rapidly acc elerates charged pa rticles, with electrons reaching 
(9(1) MeV and ions 0{1) GeV/nucleon ( lEmslie et al.ll2004| ). Explaining this energization 
requires accounting not only for the relevant energy and time scales but also the resulting 
spectra, which exhibit a common shape for almost all ion species. At the same time, high 
mass-to-charge {M/Q) ions are greatly over-represented in flares, with abundances as much 
as tw o orders of magnitude higher than normal coronal values (IMason et al.l Il994j : iMason 
2007h . 



^Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218; 
kknizhni(Q)pha.jhu.edu 

^University of Maryland, College Park, MD 20742; swisdak@umd.edu, drake@umd.edu 



-2- 



In impulsive flare models, magnetic reconnection is the ultimate energy source, and so 
it is natural to consider theories in which reconnection also plays a role in particle acceler- 
ation. Some m odels, including those that rely on interactions with mag netohydrodynamic 
(MHD) waves jMilleilllQQsi iPetrosian fc Liul 120041 ) or shock acceleration (lEUison &: Ramatyl 
19851 : ISomov fc Kosud EqOtI ) . use reconnection only as an indirect source that provides an 
environment in which other energization processes can occur. In other models, the mag- 
netic energy released during reconnect ion is more di rectly channeled to particles through 
various processes — DC e lectric fields (jHolmanlll985l ). int eractions with mu ltiple magnetic 



islands ( lOnofri et al.ll2006l ). first -order Fer mi acceleration ( jPrake et al.ll2010l ). or the pickup 
of collisionally ionized neutrals (IWulll996l ). 



Also a member of this latter group is the direct heating of ions in reconnection exhausts 
( iKrauss-Varban &: Welschl l2006l : iDrake et al.l l2009bl ). in which the perpendicular and par- 
allel (relative to the magnetic field) temperatures of ions jump after traversing the narrow 
boundary layer separating the ions that fiow slowly in from upstream from the reconnec- 
tion exhaust, which travels at the Alfven speed Bj^^i^p where B is the strength of the 
magnetic field and p is the density. However these works considered the weak guide fielcil] 
limit in which ion heating is para llel, rather than transverse, to the local magnetic field; 
Cranmer &: van BallegooiienI (120031) have shown that in the extended solar corona, Ti ^ Ty. 
Subsequently, Drake et al.l ( 2009a ) used test particle^ in a Hall MHD simulation with a large 
guide field (five times larger than the reconnecting field) to confirm that ions above a critical 
value of MjQ become demagnetized. They suggested that ions crossing into reconnection 
outflows can become non-adiabatic, and hence behave like pickup particle^ (jMobius et al. 



19851 ). while gaining an effective thermal velocity equal to the Alfven speed and deriv ed a 
Ml Q- based threshold for this behavior. This process is similar to an earlier proposal by IWu 
(]l996l ) that ion acceleration in impulsive flares can occur via reconnection-associated pickup, 
although in that case the accelerated ions wer e produced by neut ral-particle ionization in 
the lower corona. Later hybrid simulations by IWang et al.l (1200 ll ) conflrmed that injected 
protons (mimicking newly ionized particles) did behave like pickup particles in this scenario. 

In this Letter, we use a kinetic particle-in-cell (PIC) simulation to track two types of 
ions self-consistently, i.e. without resorting to test particles, to determine whether particles 



^ A guide field is a component of tlie magnetic field perpendicular to the reconnection plane. Most coronal 
reconnection is guide field reconnection. 

^Test particles are particles that move under the influence of the simulation's electromagnetic fields, but 
do have any self-consistent effect on the computation. 

•^The pickup process refers to the ionization of a neutral atom with velocity « embedded in a high- 
velocity plasma. It plays an important role throughout the heliosphere, and particularly in the solar wind. 
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above the critical value of M/Q behave like pickup particles in d ynami c electromagnetic 
fields. Ions with M/Q below the threshold derived in iDrake et al.l (l2009al ) (protons, in this 
case) are adiabatic and undergo very little heating as they move between the upstream 
plasma and the reconnection exhaust, while particles above the threshold (a particles) gain 
an effective thermal velocity equal to the exhaust velocity after crossing the narrow boundary 
layer surrounding the exhaust. 

The transition between adiabatic and non-adiabatic behavior depends on the ratio be- 
tween a particle's c yclotron period and the the time it takes to cross the boundary layer 
(IDrake et al.ll2009al ). An adiabatic particle turns sharply in the outfiowing direction upon 
entering the exhaust, conserving its magnetic moment /i = m6vj_/2B, where Sv± is the ion 
perpendicular velocity with the E x B contribution subtracted (the ion perpendicular veloc- 
ity is taken relative to B). However, particles which behave non-adiabatically move in the 
direction of the local electric field upon entering the exhaust and not in the direction of the 
local E X B velocity. The sudden change from slow upstream infiow to downstream Alfvenic 
outflow causes particles with high M/Q to see a jump in their magnetic moments. 



2. NUMERICAL SIMULATIONS 



We carry out simulations using the code p3d (IZeiler et al.ll2002l ). Like all PIC codes, it 
tracks individual particles (~ 10^ in this work) as they move through electromagnetic fields 
that are defined on a mesh. Unlike more traditional fiuid representations (e.g., MHD), PIC 
codes correctly treat small lengthscales and fast timescales, which are particularly important 
for understanding the x-line and separatrices during magnetic reconnection. 

The simulated system is periodic in the x — y plane, where fiow into and away from the 
x-line are parallel to y and x, respectively, and the guide magnetic field and reconnection 
electric field parallel z. Th e initial magnetic field and density profiles are based on the Harris 
equilibrium ( HarrisI 1962 ). The reconnecting magnetic field is given by = tanh[(?/ — 



Ly/4:)/wo] — tanh[(|/ — 3Ly/4:)/wo] — 1, where Wq and Ly are the half- width of the current 
sheets and the box size in the y direction. The density comprises an ambient background 
and two current sheets in which the density rises in order to maintain pressure balance with 
the magnetic field. We initiate reconnection with a small initial magnetic perturbation that 
produces a single magnetic island on each current layer. 

The code is written in normalized units in which magnetic fields are scaled to the 
asymptotic value of the reversed field Bq^, densities to the value at the center of the re- 
connecting current sheet minus the uniform background density, velocities to the proton 
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Alfven speed ca = -Box/ A/47rmJ^, times to the inverse proton cyclotron frequency in Bq^.-, 
VL~^ = nipc/eBox, lengths to the proton inertial length dp = ca/^px and temperatures to 



2 



The proton to electron mass ratio is taken to be 25, in order to minimize the difference 
between pertinent l ength scales and hence run as larg e a simulation do main as possible. 



It has been shown ( jShay et al.l Il998l : iHesse et al.l Il999l : IShay et al.l 120071 ) that the rate of 



magnetic reconnection and structure of the outflow exhaust do not depend on this ratio, and 
neither, therefore, does the ion heating examined here, which depends only on the exhaust 
geometry. The simulation assumes d/dz = 0, i.e. that field and particle quantities do not 
vary in the out-of-plane direction, making this a two-dimensional simulation. 

In addition to the usual protons and electrons, we also include a number density of 1% 
^He"'""'" (a) particles in the background particle population and gave them an initial temper- 
ature equal to that of the protons. This number density does not affect the reconnection 
dynamics appreciably, while still providing a large sample of particles with M/Q > 1, where 
M/Q is normalized to the proton value. Each particle (protons and a's) is assigned a unique 
tag number, allowing individual particles to be tracked throughout the simulation. 

In Fig. [1] we show an overview of results from a simulation with a computational domain 
Lx X Ly = 102.4 X 51.2 dp and an initial guide field Bq^ = ^.OBqx at t = 200^2"^^. The grid 
spacing for this run is 0.025 dp, the electron, proton, and a temperatures, Te = Tp = Ta = 
0.25'mpc\, are initially uniform and the velocity of light is ISc^- The half-width of the initial 
current sheet, Wq, is 1 dp and the background density is 0.2. Panel (a) depicts the total 
out-of-plane current density Jz centered around the x-line (at x/dp ^ 32 and y/dp ~ 13) of 
one of the current sheets. Magnetic field lines (not shown) roughly trace contours of Jz. 

Ambient plasma from above and below the current sheet slowly flows toward the current 
sheet while embedded in oppositely directed magnetic field (to the right above the layer, to 
the left below). Reconnected field lines are highly bent and, to reduce their magnetic tension, 
rapidly move away from the x-line, dragging plasma with them. Panels (b) and (c) show 
the proton and a outflow velocities Vp^ and Vax- The similarity between the two makes 
it clear that both the protons and the a's participate in the reconnection outflow which, 
outside of the immediate vicinity of the X-line, has a magnitude of Vx ~ cEy/Bz ~ ca 
a 1, in our normalized units). A comparison of Fig. [1] to frames (a) and (b) of Fig. 1 in 



Drake et al.l (j2009al ) (which shows results from a run otherwise identical but for the presence 
of the a particles) demonstrates that the a's do not significantly change the structure of the 
reconnection exhaust. 
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3. ION PICKUP AND HEATING 

Particle acceleration is controlled by the structure and magnitude of the electric field. 
During reconnection a strong transverse electric field Ey = —EzBqz/ By develops in the 
exhaust to force E ■ B = 0; its structure to the left of the x-line is shown in the background 
of Fig. [2](a). Particles enter the exhaust with a velocity Vy ~ O.lc^. Any energy gain 
is determined by whether particles crossing the exhaust boundary, which has scale length 
given by the ion sound Larmor radius ps (= Vs/^lp, where Vg is the plasma sound speed), are 
adiabatic. Non-adiabatic particles cross the boundary in a time Tc that is short compared 
with their cyclotron period, ~ Ps/vy ~ lOp^/c^ < n/Qpx, or 

(1) 




where f3px = StttiT/Bq^ ( iDrake et al.ll2009al ). Thus, in the present simulations, where the 
upstream (3px = 0.2, equation [1] gives M/Q > 1 for non-adiabatic behavior and so protons are 
marginally adiabatic while a's {M/Q = 2) are not. Since -Ey < in this case, non-adiabatic 
positively charged ions entering the exhaust from below will be pushed out immediately, 
preventing them from being caught up in the exhaust. However, non-adiabatic positively 
charged ions entering from the top will find themselves essentially at rest in the simulation 
frame while the outflow moves past at roughly the Alfven speed. Such particles will undergo 
an E X B drift, but with a "thermal velocity" equal to the Alfven speed and have trajectories 
resembling cycloids. This process is analogous to that undergone by stationary neutral atoms 
surrounded by the moving solar wind. If ionized, the new ion first moves in the direction of 
the motional electric field in order to gain the necessary energy to fiow with the rest of the 
wind. As it gets "p icked up", it gains a thermal velocity equal to the solar wind velocity 



(iMobhis et al.lll985l 



We randomly selected 500 protons and 500 a particles from the 7.5 dpX3dp box upstream 
of the exhaust shown in Fig. [2](a) at t = 200Q~^ and followed their trajectories for 25 
In Fig. [2]^a) we plot a representative trajectory for a proton, shown in black, and an a 
shown in green, over a background of Ey. (Note that the overlaid trajectories in (a) are 
calculated in the fully self-consistent simulation, while the background of Ey is a snapshot 
from t = 202^2"^?^.) The proton, which remains adiabatic, immediately moves downstream 
upon entering the exhaust, while the a particle moves in the direction of Ey before being 
picked up by the E x B drift. Panel (b) displays the time evolution of the proton (black) 
and a (green) magnetic moments (scaled by mass). The vertical red line corresponds to the 
time at which Ey is shown in (a). After crossing the boundary layer into the exhaust, the 
a becomes demagnetized, as indicated by the jump in p, a trend seen for all of the tracked 



- 6 - 



a's. In (c) we plot the magnetic moments of all 500 protons (black) and all 500 a's (green) 
after entering the exhaust versus their moments at t = 200Q~^. For each particle /ifinai was 
measured when the particle crossed a specified horizontal position at the downstream edge 
of the exhaust, around 6 dp in Fig. I2^a). For reference, we overplot a line of unit slope, 
which corresponds to exact /i conservation. The clustering of protons near this line and 
large values of /i/m reached by the a's clearly shows the adiabatic nature of the former and 
the non-adiabatic nature of the latter. 

In Fig. [3] we show the temperatures of both species. Panels (a) and (b) depict the 
perpendicular (to the magnetic field) a and proton temperatures, while (c) and (d) show the 
parallel temperatures. The a temperature increase is greater than the proton temperature in 
the perpendicular direction (note the different color bar scales in the two panels). Indeed, the 
temperature increase of the a's is na ore than mass proportional, consistent with observations 



(jCranmer fc van Ballegooijenll2003l ). This is also evident in frames (e) and (f), which are 
cuts through the perpendicular and parallel temperature plots for the a's (red) and protons 
(black). The weak heat ing of the protons is consistent with the adiabatic behavior shown in 
Fig. 121 The analysis of iDrake et al.l (l2009al ) predicts that, with a guide field of Bq^ = 2Bqx, 



the proton temperature will change by 



^^11 = 1?^' ~ T' = (2) 

in the exhaust, and the a temperature will change by: 

AT,| = 0; ATx = ^m^vl ~ 2vl. (3) 

For f ^ ~ 2 (see Fig.[T]) these jumps are in reasonable agreement with the observed variations. 
Differences from the predicted values, in particular the changes in T± for the protons and Ty 
for the a's, presumably arise from corrections to equations ([2]) and ([3]) due to a mixture of 
adiabatic and non-adiabatic behavior by the particles. 

In Fig. m we show how the velocity distributions of the a's and protons change as they 
move from the upstream to the downstream region. Panels (a) and (d) depict the upstream 
a and proton velocity distribution in the Vx — Vy plane. The protons and a's were given the 
same initial temperature, so upstream of the exhaust, the protons' mean thermal velocity 
is higher than the a particles'. The small negative Vy component upstream in both species 
shows the inflow toward the reconnection exhaust. After crossing the narrow boundary layer, 
the a's get picked up by the Alfvenic outflow (at this time the local density is ~ 0.2, so the 
local Alfven speed is ~ "\/5), and their thermal velocity increases much more than that of the 
adiabatic protons. The downstream velocity distributions for the a's and protons are shown 
in the v^—Vy plane in (b) and (e) and in the v^ — Vz plane in (c) and (f), both calculated inside 
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a box located between 5 — 12.5 dp in the x-direction and 11.4 — 13.4 dp in the y-direction. 
Since the dominant B component is the guide field, and Vy are essentially perpendicular 
velocity components, and the parallel velocity. The protons exhibit very little heating in 
the Vx — Vy plane (Fig.|3](e)), consistent with adiabatic behavior. They are modestly heated in 
the z-direction (Fig. 111(f)), consistent with Ty in Eq. [2l The a particles are strongly heated in 
the y-direction (Fig. |4](b)) and are beginning to form a ring distribution that is characteristic 
of pickup behavior. There is modest heating of a particles in the z-direction (Fig. Hl^c)), but 
the similar structure in (c) and (f) suggests that the heating mechanism is the same for both 
protons and a's, and it is possible that the a particles are not completely non-adiabatic. 



4. DISCUSSION 



Using self-consistent trackin g of particle traject ories, we have shown that ions above a 
critical mass-to-charge threshold (jPrake et al.ll2009al ) behave like pickup ions in reconnection 
exhausts, with /i changing due to a sharp increase in Vj_. Ion energy increments of ~ 25 
keV/nucleon are predicted for typical coronal parameters of -B = 50 G and n = 10^/cm'^. 
Ions below this threshold are only weakly heated. This transition only exists for reconnection 
with a guide field which, however, is the typical case in the solar corona. Coronal observations 
have revealed that the abundances of high mass-to-charge ions are enhanced in solar flares, 
with the strength of the enhancement depending only on M/Q. The fact that we observe 
non-adiabatic behavior and associated strong heating for particles with M/Q > 1, while 
the proton heating remains weak, suggests that reconnection might explain the abundance 
enhancements in impulsive flares. Abundance enhancements should occur because high M/ Q 
ions are heated at lower values of the reconnecting magnetic fleld strength (see Eq. [T]) than 
protons. Furthermore, the increase in T\ /T\\ in the exhaust seen here is co nsistent with 
that observed in the extended solar corona (jCranmer fc van Ballegooijenll2003l ). although it 
should be noted that the number density of a particles used here is slightly less than what 



is observed in the corona 



n. 



a, corona 



5% - 10%). 



Observatio ns near 1 AU of sol ar wind reconnection events with the Wind spacecraft (as, 
for example, in iPhan et al.l ( 120101 )) should be able to measure Ty and T± for both protons 
and a particles in order to test the mechanism suggested in this work. Moreover, provided 
its instrumentation can differentiate between different M/Q ions, the upcoming Solar Probe 
Plus mission, with a planned perihelion of ~ 9Rq (which lies within the outer corona), 
should also provide an excellent test of these predictions. 



Finally, it is widely believed that some process converts a fraction of the energy found 
in the convective motions of the solar photosphere into the heat that ensures the continuous 
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existence of a O{10^ K) corona and accelerates the solar wind. Broadly speaking the two 
most likely candidates are wave heating — in which oscillations generated in the photosphere 
travel into the corona, develop into turbulence, and eventually dissipate — and magnetic 
reconnection, in which the topological reorganization of the magnetic field releases energy 
and heats the plasma. 

Measurements by the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) 
and Ultraviolet Coronal Spectrometer (UVCS) instruments of the SOHO (Solar and He- 
liospheric Observatory) spacecraft provide significant constraints on any theory of coronal 
heating. In particular, at heights of 2 — 3Rq protons have a slight temperature anisotropy 
(in the Tj_ > Tn sense) while heavier ions (represe nted by 0^~^) are strongly anisotropic, 



with T^/Tn > 10 (jCranmer &: van Ballegooijen 



20031 ). Interestingly, the process discussed in 



this work should be active in the region in question and produces temperature anisotropics 
consistent with these results. 
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Fig. 1. — Overview of a PIC simulation with an initial guide field Bqz = 2Box- Panel (a): 
the total out-of-plane current density J^; panel (b): the proton outflow velocity Vp^; panel 
(c): the a particle outflow velocity Vax 
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Fig. 2. — Panel (a): Trajectories for a proton (black) and a particle (green) randomly picked 
from a 7.5 x 3cip box (shown in white) are overlaid on a snapshot of Ey. Panel (b): The 
magnetic moment per mass as a function of time for the two particles in (a). The red line 
represents the time of the snapshot of Ey. Panel (c): For 500 protons (black) and 500 
a's (green) selected at random from the white box in (a), their final magnetic moments 
plotted against their initial magnetic moments. The red line has unit slope and represents 
the expected result if /i/m were invariant. 



- 12 - 




(e) 10 20 30 40 50 




0.0 



5 10 15 20 25 

y/dp 

Fig. 3. — Spatially smoothed temperature components at t = 200ri~^^. Panel (a): T± for 
a particles; panel (b): Tj_ for protons; panel (c): Ty for a's; panel (d): Ty for protons. 
Panels (e) and (f): Cuts through the exhaust of the perpendicular and parallel components, 
respectively, for a's (red) and protons (black). The locations of the cuts are shown by the 
vertical lines in (a)-(d). 
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Fig. 4. — 2D velocity distribution functions upstream and downstream of the exhaust. Panels 
(a) and (d): upstream Vx — Vy distributions for a's and protons, respectively. Panels (b) and 
(e): downstream — Vy distributions for a's and protons; panels (c) and (f): downstream 
Vx — Vz distributions. 



